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We search for standard model single-top-quark production in the VF-gluon fusion and W* channels 
using 106 pb" 1 of data from pp collisions at = 1.8 TeV collected with the Collider Detector at 
Fermilab. We set an upper limit at 95% C.L. on the combined W-gluon fusion and W* single-top 
cross section of 14 pb, roughly six times larger than the standard model prediction. Separate 95% 
C.L. upper limits in the W-gluon fusion and W* channels are also determined and are found to be 
13 and 18 pb, respectively. 

PACS numbers: 14.65.Ha, 12.15.Ji, 13.85.Rm 



The observation of the top quark in pp collisions at the 
Fermilab Tevatron has relied on pair production through 
the strong interaction, typically qq — > ti. A top quark 
can also be produced singly, in association with a b quark, 
through the electroweak interaction jl[]. The two domi- 
nant "single-top" processes are 'Wg"(i.e. W-gluon fu- 
sion, qg — > tbq') and U W*" (qq' — > tb). Within the con- 
text of the standard model, a measurement of the rate of 
these processes at a hadron collider allows a determina- 
tion of the Cabibbo-Kobayashi-Maskawa matrix element 
V t b §■ Assuming |V t fc| = 1, the predicted cross sections 
for Wg and W* are 1.7 pb Q and 0.7 pb || respec- 
tively, compared to 5.1 pb for tt pair production The 
D0 Collaboration has recently published 95% confidence 
level (C.L.) upper limits of 22 pb on Wg and 17 pb on W* 
production In this Letter we report on two searches, 
one for the two single-top processes combined, and the 
other for each process separately. 

The expected final state of a single-top event consists of 
W-decay products plus two or more jets, including one b- 
quark jet from the decay of the top quark. In W* events, 
we expect a second 6-quark jet from the W*tb vertex. In 
Wg events, we expect a jet originating from the recoiling 
light quark and a second 6-quark jet produced through 
the splitting of the initial-state gluon. This 6-quark jet 
is produced at larger absolute value of pseudorapidity 
and lower transverse momentum than in W* events 

Single-top processes are harder to observe than ti pro- 
duction because their cross section is smaller and their 
final state, containing fewer jets, competes with a larger 
147+multijet background from QCD. A priori we do not 
expect sensitivity to the standard model cross section in 
the presently available data. However, a number of new 
physics processes could enhance the single-top produc- 
tion rate, motivating a search ||, 0]. 

Our measurement uses 106 ± 4 pb" 1 of data from pp 
collisions at y/s = 1.8 TeV collected with the Collider 
Detector at Fermilab between 1992 and 1995 ("Run I"). 
The detector is described in detail elsewhere We 
restrict our single-top search to events with evidence of a 
leptonic VF-decay: an isolated |ll[] electron (muon) can- 
didate with E T (P T ) > 20 GeV (GeV/c) and missing 



transverse energy [[L2j fir > 20 GeV from the neutrino. 
We remove events that were identified in a previous CDF 
analysis |L3| as ti dilepton candidates. Events with a sec- 
ond, same-flavor and opposite-charge lepton that forms 
an invariant mass with the first lepton between 75 and 
105 GeV/c 2 are rejected as likely to have come from Z° 
boson decays. Furthermore, to reject those dilepton ti 
or Z° candidates where one lepton fails our electron or 
muon identification, we also remove events that contain 
a track with Pt > 15 GeV/c and charge opposite that 
of the primary lepton, and such that the total Pt of all 
tracks in a cone of radius AR = \J A77 2 + A</> 2 = 0.4 
around this track is less than 2 GeV/c Q|. Jets are 
formed as clusters of calorimeter towers within cones of 
fixed radius AR = 0.4. Events are required to have one, 
two, or three jets with Et > 15 GeV and \rj\ < 2.0; 
at least one jet must be identified as likely to contain a 
b quark ("6-tagged") using displaced- vertex information 
from the silicon vertex detector (SVX) 14 . If a second 



jet in the event is also 6-tagged, either in the SVX or by 
the presence of a soft lepton indicative of semilcptonic 
b decay, the event is labeled "double-tag" , otherwise it 
is labeled "single-tag". The above event selection cuts 
are common to our combined and separate searches for 
the two single-top processes. Additional cuts are applied 
within each analysis. 

We first describe our search for single-top production 
in the Wg and W* channels combined. The expected 
signal significance is improved by requiring the invariant 
mass Mfvb reconstructed from the lepton, neutrino, and 
highest-_ET ^-tagged jet, to lie in a window around the top 
quark mass, 140 < M (l/b < 210 GeV/c 2 . The neutrino 
momentum is obtained from the fir and the constraint 
that Mg v = Myy |lq |. The variable discriminates 
against both non-top and ti backgrounds, in the latter 
case because combinatorial errors in assigning partons to 
final-state jets broaden the distribution compared 

to single top. 

We determine the efficiency of our selection criteria 
from events generated by the PYTHIA Monte Carlo pro- 
gram Q and subjected to a CDF detector simulation. 
The acceptance times branching ratio is (1.7 ± 0.3)% for 
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Combined Search 


Separate 


Search 




W + 1,2, 3 jets 


W + 2 jets 


Process 




Single-tag 


Double-tag 


Wg 


3.0 ±0.6 


1.4 ±0.3 


0.04 ±0.01 


W* 


1.3 ± 0.2 


0.55 ±0.15 


0.32 ±0.06 


tt 


8.4 ± 2.7 


1.4 ±0.5 


0.7 ±0.2 


non-top 


54 ± 12 


10 ± 2 


1.6 ±0.4 


Total 


67 ± 12 


14 ± 2 


2.7 ±0.5 


Observed 


65 


15 


6 



TABLE I: Expected numbers of signal and background events 
passing all cuts in the W+jets data sample, compared with 
observations. The uncertainties on the expected numbers of 
single-top events do not include uncertainties on the theoret- 
ical cross section calculations. 



each of the two single-top processes. The largest contri- 
butions to the acceptance uncertainties come from lepton 
triggering and identification (10%), and ^-tagging (10%). 
Combining these acceptances with the cross sections pre- 
dicted by theory j| [| and the size of the CDF Run I 
dataset, we expect a total signal yield of 4.3 events. 

Expectations for signal and background rates are listed 
in the second column of Table |. We estimate the tt 
background from a herwig Monte Carlo calculation ]Tij ] 
followed by a detector simulation. Normalizing to the 
theoretically-predicted cross section, er t j = 5.1 ± 0.9 pb 
||, we expect 8.4 ± 2.7 events to survive our selection 
criteria, where the uncertainty includes theoretical and 
acceptance contributions. 

The largest component of the non-tt background in the 
SVX-tagged VF±jets sample is inclusive W production 
in association with heavy-flavor jets (e.g. pp — » Wg, fol- 
lowed by g — > bb). Additional sources include "mistags," 
in which a light-quark jet is erroneously identified as 
heavy flavor, "non-IF" (e.g. direct bb production), and 
smaller contributions from WW, WZ, and Z±heavy- 
flavor Jl4j]. The mistag and non-IF rates are estimated 
from data, the W+heavy-flavor rates from Monte Carlo 
normalized to data, and the smaller sources such as dibo- 
son production from Monte Carlo normalized to theory 
predictions Q . The total non-top background expecta- 
tion is 54±12 events. The uncertainty on our background 
includes the effect of varying the top mass by its uncer- 
tainty of ±5 GeV/c 2 . 

To measure the combined Wg + W* single-top produc- 
tion cross section, we use a kinematic variable whose dis- 
tribution is very similar for the two single-top processes 
and is different for background processes: the scalar sum 
Ht of Ipr and the transverse energies of the lepton and 
all jets in the event. We perform an unbinned maximum- 
likelihood fit of the Ht distribution from data to a lin- 
ear superposition of the expected Ht distributions from 
single-top signal, tt and non-top backgrounds. We model 
the shape of the Ht distribution for all sources of non-top 
background with VECBOS-generated |l8| events contain- 




100 150 200 250 300 350 400 450 500 
H T (GeV) 

FIG. 1: The Ht distribution for data in the combined search, 
compared with smoothed Monte Carlo predictions for signal 
and backgrounds (second column in Table |) . Ht is the scalar 
sum of Ipr and the transverse energies of the lepton and all jets 
in the event. The inset shows that the Monte Carlo modeling 
of Ht is very similar for both signal processes. 



ing a W plus two partons that we force to be a bb pair. We 
have checked that VECBOS reproduces the Ht and Mg^ 
distributions for the 6-tagged W + 1-jet data before the 
Mgvi, cut, a sample in which the non-top backgrounds are 
expected to dominate. In the search sample, the observed 
Ht distribution agrees with the spectrum derived from 
Monte Carlo calculations when the latter are normalized 
to the a priori predicted numbers of events (Figure [l]) . 

We set an upper limit on the cross section using the 
likelihood function: 

C(P s ,P t t,Pnt) = Gi(/3 tt -) X G 2 (Pnt) X C s ha,pe(Ps,PtuPnt), 

where (3 S , (3 t i and (3 n t are fit parameters representing, 
respectively, factors by which the standard model cross 
section predictions for single-top, tt and non-top must 
be multiplied to fit the data. The functions Gi and G2 
are Gaussian densities constraining the background fac- 
tors P t t and P n t to unity, and £ s hapc represents the joint 
probability density for observing the A^bs data events at 
their respective values of Ht'- 



N obs - 

Met e 



N ohs \ 



N 



x tj PsF s (Hn) + f3uF t - t {H Tl ) + (3 n tF nt (H T i 



Mfit 
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Wg + W* 


Wg 


W* 




Source 


Si i 


AS 


fin 


AS 


fin 


AS 




TpT Rt QPfllp 


0.01 


0.25 


0.01 


0.02 


0.01 


0.06 




Initial-state radiation 


0.02 


0.15 


0.06 


0.07 


0.06 


0.13 




Final-state radiation 


0.03 


0.02 


0.07 


0.02 


0.05 


0.01 




Part on distributions 


0.04 


0.02 


0.01 


0.03 


0.01 


0.02 


o 


Signal generator 


0.02 


0.25 


0.08 


0.03 


0.07 


0.12 




Background model 




0.04 




0.12 




0.18 


a -1 


Top mass 


0.04 


0.01 


0.01 


0.12 


0.00 


0.35 


c 


Trigger & lepton id. 


0.10 




0.10 




0.10 




> 

Ld 


6-tag efficiency 


0.10 




0.10 




0.10 






Luminosity 


0.04 




0.04 




0.04 






Total 


0.16 


0.39 


0.19 


0.19 


0.18 


0.44 





TABLE II: Systematic uncertainties on the fit result for j3 s 
in the combined search (Wg + W*), and for f3wg and fiw* in 
the separate search (see text). The 5n columns list fractional 
uncertainties due to signal normalization effects and the AS 
columns absolute uncertainties due to effects on the shapes of 
the fitted distributions. 



In this expression, fiR t = ftsfJ-s + PtifJ-ti + PntUnt, where 
/i s , jtt t f and [i n t are the predicted numbers of events, 
and F{Ht) are smoothed Ht distributions for signal and 
background, normalized to unity. The maximum of C is 
obtained for j3 s = 2.0 ± 1.8, where the uncertainty is sta- 
tistical only and includes the effect of correlations with 
the other fit parameters. 

To extract Bayesian upper limits on the single-top 
production rate, we construct a probability distribution 
fifis) by maximizing C((3 S , /%, (3 n t) with respect to /3 t j 
and f3 n t for each value of (3 S , and multiplying the result 
with a flat prior distribution for [3 S . We then convolute 
f(f3 s ) with two Gaussian smearing functions. The first 
one has width (3 s 5n, where Sn is the sum in quadrature 
of all the normalization uncertainties listed in Table |H|. 
The width of the second smearing Gaussian is the sum 
in quadrature of all the systematic uncertainties relative 
to the shape of the Ht distribution (AS in Table |l|) . Fi- 
nally, the smeared distribution is integrated to find the 
95% C.L. upper limit on single-top production. We find 
this limit to be /3 S 95 = 5.9, corresponding to a cross sec- 
tion of 14 pb. 

Because of significant differences in the final-state kine- 
matics of the two single-top processes, it is possible to 
search for them separately. This is interesting, because 
an exotic single-top production mechanism may con- 
tribute to one and not the other, for example a heavy W' 
decaying to a tb quark pair adding to the apparent W* 
rate For the separated search, we use events in the 
TF+2-jets sample only and consider two non-overlapping 
subsamples. The first one consists of single-tag events 
in which the reconstructed top mass lies in the window 
145 < Mt v b < 205 GeV/c 2 , and the second consists of 
double-tag events. The expected compositions, calcu- 
lated in the same way as for the combined analysis, are 
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FIG. 2: Top: distribution of the product Q x r\ of the lep- 
ton charge and the untagged jet pseudorapidity for single-tag 
W+2-jets events. Bottom: distribution of the reconstructed 
top mass for double-tag events. The data are compared with 
expectations for signal and backgrounds (third and fourth 
columns in Table ffl). 



shown in the last two columns of Table |: in the single- 
tag sample, Wg is about 2.5 times larger than W*\ in 
the double-tag sample, W* is about 7.5 times larger than 
Wg. 

The Wg component in the single-tag sample can be 
measured by considering that the light-quark jet in 
Wg events is about twice as likely to be in the same 
hemisphere as the outgoing (anti)proton beam when a 
(anti)top quark is produced. Thus the product Q x 77 of 
the primary lepton charge and the untagged jet pseudo- 
rapidity has a strongly asymmetric distribution. In the 
double-tag sample, the W* component can be extracted 
from the distribution of M^. In this case, since both jets 
are tagged, the 6-jet with the largest rj (—77) is used in 
forming for a t (t) decay, as determined by the sign 
of the primary lepton in the event, an assignment that is 
expected to be correct 64% of the time. The Q x r\ and 
Mivb distributions for the data are compared to expec- 
tations for signal and background in Figure |^. For the 
separate Wg and W* searches, we use a herwig Monte 
Carlo calculation to model our signals. 

A binned maximum-likelihood fit is used to extract the 
amounts of Wg and W* present in the T / F+2-jets data. 
The likelihood function has the following form: 

Gi(Atl) X G 2 (/?„ t i) X &(j3wg, 0W* , foil, Pntl) X 



6 



where the fit parameters are factors by which the pre- 
dicted numbers of Wg {(3wg), W* (@w*), single-tag ti 
(/Stti), double-tag ti (/3tt2), single-tag non-top ((3 n ti) and 
double-tag non-top {j3 n n) events must be multiplied to 
fit the data. The Gi functions are Gaussian constraints 
on the normalizations of the various backgrounds, C\ is 
a binned Poisson likelihood for the Q x 77 distribution of 
single-tag events, and £2 is a binned Poisson likelihood 
for the distribution of double-tag events. 

The result of the maximum- likelihood fit for the single- 
top content of the data is — 0.61 4 'q Wg events and 7.6^|' g 
W* events. The systematic uncertainties are listed in Ta- 
bic [n| We extract upper limits on the individual single- 
top processes in the same way as for the combined search. 
At the 95% C.L., we find upper limits of 13 and 18 pb 
on single-top production in the Wg and W* channels, 
respectively. These two limits are correlated since they 
are derived from the same likelihood function. 

In summary, we conclude that electroweak tb produc- 
tion is out of reach in the Run I CDF data set. At 
the 95% C.L., we set an upper limit on the combined 
Wg + W* single-top cross section of 14 pb. Separate 
95% C.L. upper limits in the Wg and W* channels are 
13 and 18 pb, respectively. 
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